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If magnetic semiconductors are ever to find 
wide application in real spintronic devices, their 
magnetic and electronic properties will require 
tailoring in much the same way that band gaps 
are engineered in conventional semiconductors. 
Unfortunately, no systematic understanding yet 
exists of how, or even whether, properties such 
as Curie temperatures and band gaps are related 
in magnetic semiconductors. Here we explore 
theoretically these and other relationships within 
64 members of a single materials class, the Mn- 
doped II-IV-V2 chalcopyrites, three of which are 
already known experimentally to be ferromag- 
netic semiconductors. Our first-principles results 
reveal a variation of magnetic properties across 
different materials that cannot be explained by 
either of the two dominant models of ferromag- 
netism in semiconductors. Based on our results 
for structural, electronic, and magnetic proper- 
ties, we identify a small number of new stable 
chalcopyrites with excellent prospects for ferro- 
magnetism. 

One of the tantalizing promises of dilute magnetic 
semiconductors is to conjoin, in a single material, the 
advantages of nonvolatility and bandgap engineering^'^'"^. 
In normal (nonmagnetic) semiconductors, precise tailor- 
ing of band structure is possible due to our detailed 
understanding — both empirical and theoretical — of the 
mechanisms underlying band formation. Unlike the sit- 
uation for normal semiconductors, however, there does 
not yet exist any comprehensive understanding of how 
the magnetic properties of magnetic semiconductors are 
related to their structural and electronic properties. 

For example, there is no firm understanding of the rela- 
tionship between the properties of the host semiconduc- 
tor and the Curie temperatures attainable by doping with 
magnetic impurities such as Mn. Experimentally, a num- 
ber of materials issues — including Mn content, compen- 
sation, and phase purity — are not yet sufficiently under 
control to permit a systematic description across differ- 
ent semiconductors. Theoretically, there is an emerging 
consensus that while the ferromagnetism originates from 
the interaction between itinerant electrons (or holes) and 
localized electrons, there is no single model that can sys- 
tematically relate the physical and electronic properties 
of the host semiconductor to its resulting magnetic prop- 
erties when doped^. Models based on the weak- and 
strong-coupling limits of this interaction (the Zener'^'^ 
and double-exchange^'^ models, respectively) do make 
specific predictions, but it is increasingly evident that 
real materials are not sufficiently well described by either 



limit for these predictions to be reliable across different 
hosts. 

Here we address this problem by exploring theoreti- 
cally these relationships within a single materials class: 
the II-IV-V2 chalcopyrites consisting of the 64 possible 
combinations of II=(Be,Mg,Zn,Cd), IV=(C,Si,Ge,Sn), 
and V=(N,P,As,Sb). Three of these— CdGePa, ZnGePz, 
and ZnSnAs2 — were recently shown to become ferro- 
magnetic upon Mn doping, with remarkably high Curie 
temperaturesSiifliiiiiS. It is an open question whether 
other chalcopyrites can be made ferromagnetic; nor is it 
known which properties of the host chalcopyrite are im- 
portant for ferromagnetism and, in particular, for high 
Curie temperatures. These are the central issues we ad- 
dress in this work. 

For each chalcopyrite, we use density-functional theory 
(DFT) in the generalized-gradient approximationiSii^ii^ 
to compute several properties of the undoped hosts: equi- 
librium crystal structure, i.e. lattice parameters a and c 
and the single internal coordinate; phase stability, as de- 
termined by the enthalpy of formation; and electronic 
band gap, which may be direct or indirect. A compre- 
hensive survey of the magnetic properties within such 
a large class of materials is impractical, so we consider 
here only a few key quantities that control magnetism. 
Since Mn is isoelectronic with the group-II elements, it 
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FIG. 1: Theoretical band gap vs. lattice constant for the II- 
IV- V2 chalcopyrites considered here. Shaded areas indicate 
hosts that are expected to be closely (±2%) lattice matched 
to either GaN, Si, GaAs, or InAs. Filled symbols denote sta- 
ble host compounds (negative enthalpies of formation) . Chal- 
copyrites whose calculated gaps are zero (see text) are not 
shown. 
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is believed that carrier-mediated ferromagnetism in the 
II-IV-V2 chalcopyrites is favored when Mn occupies the 
group-IV site, where it is hkely to be electrically activei^. 
For this reason we first explore the competition between 
substitutional Mn doping on the II and IV sites by com- 
puting the Mn impurity formation energy on each site. 
Second, we address whether occupation of the group-IV 
sites is indeed likely to lead to a ferromagnetic phase; 
in particular, for each host we compute the Heisenberg 
spin coupling for an experimentally reasonable Mn con- 
centration. The magnitude (and sign) of this coupling 
should serve as a useful roadmap for the experimental 
exploration of new ferromagnetic chalcopyrites. 

Of the 64 chalcopyrites considered here, about 1/3 are 
known to exist and have been studied experimentally. A 
few crystallize in both the chalcopyrite and sphalerite 
structure (e.g. CdGeP2 and ZnSnSb2), while two take 
the wurzite structure (BeSiN2 and ZnGeN2); for simplic- 
ity we treat all 64 in the chalcopyrite structure. For those 
with measured lattice parameters, we can compare our 
theoretical predictions. As expected from DFT calcula- 
tions, a is given very accurately (within ~1%), while c 
is less so (^15%). As Fig. 1 shows, nearly half of the 64 
chalcopyrites enjoy a close lattice match to an important 
mainstream semiconductor. 

As with all semiconductors, DFT in the generalized- 
gradient approximation significantly underestimates 
band gaps. For chalcopyrites whose measured gaps are 
in the range 1-2 eV, the DFT values are approximately 
1 eV smaller than measured, consistent with previous 
findingsii. Hence the theoretical gaps shown in Fig. 1 
must be viewed with caution: the true gaps are likely 
to be at least 1 eV larger than predicted. Likewise, 
hosts that have no gap within DFT may in reality have 
a gap (which is likely to be 1 eV or smaller). These 
caveats aside, the trend in Fig. 1 is standard for cova- 
lently bonded semiconductors: smaller lattice constants 
are associated with larger band gaps. In particular, the 
nitrides have the smallest lattice constants and largest 
gaps, while the antimonides have the largest lattice con- 
stants and smallest gaps. From the perspective of our 
survey, the range and distribution of lattice constants and 
band gaps is satisfyingly large, guaranteeing our roadmap 
good coverage of this unexplored territory. 

Not all of the chalcopyrites considered here are ther- 
modynamically stable, even in the absence of Mn. 
Fig. 2 shows the calculated enthalpies of formation, 
Ai?/(II-IV-V2) ^ E{ll-lV-V2) - E{IT) - E{W) - 2E{Y), 
where the latter terms refer to the ground-state elemen- 
tal phases. Of the 64 chalcopyrites, 3/4 have negative en- 
thalpies, suggesting many potentially stable host materi- 
als. (We have not considered phase stability with respect 
to disproportionation into compound products, and so 
the terms "stable" and "unstable" should be understood 
here with this caveat.) The trend in Fig. 2 is similar, 
and indeed intimately related, to that of Fig. 1: smaller 
lattice constants are usually associated with higher sta- 
bility. Hence, as expected, the nitrides tend to be the 



most stable chalcopyrites and the antimonides the least 
stable. Note that there are important deviations from 
this trend: for example, nearly all the carbides are un- 
stable. 

The solubility of Mn in a host chalcopyrite is an ex- 
tremely important issue for two reasons. The first per- 
tains to the well-known tendency, within Mn-doped DMS 
materials, toward phase separation with increasing Mn 
content^'^. This tendency is reflected in the very low 
equilibrium solubility limit of Mn in, for example, IIT 
V hosts — typically of order parts per thousand or less. 
These limits can be greatly exceeded, by 2-3 orders of 
magnitude, by relying on kinetic barriers to maintain 
the metastable phase. It is not feasible to address this 
issue theoretically for the chalcopyrites. Instead, we pro- 
pose using the theoretical solubility limit of Mn as a 
rough proxy for the metastability doping limit. The 
solubility limit of Mn, Xsoi, is determined by its impu- 
rity formation energy, Ai?^(Mn), simply according to 
Xsoi — exp(— AiJ/(Mn)/fcT). In this formulation, a neg- 
ative formation energy implies that spontaneous incorpo- 
ration of Mn impurities will be limited only by kinetics^. 

The second reason concerns ferromagnetism per se. In 
general, the solubility limits for Mn substituting on the 
group-II and -IV sites will be different. Mahadevan and 
Zunger have shown that, in the case of CdGeP2, sub- 
stitutional Mnn leads to antiferromagnetic interactions 
while Mniv leads to ferromagnetic interactionsi^. Hence 
the II-IV solubility difference will greatly influence the 
stability, or even existence, of the ferromagnetic phase, 
as well as the attainable Curie temperatures. There are 
two other important issues, which we do not address here. 
One is the possibility of Mn occupying interstitial sites 
(which is increasingly likely for hosts with small lattice 
constantaiS.) or forming complexes with other Mn or na- 
tive defects; if electrically active, these can have conse- 
quences for ferromagnetism. The other is the possibility 
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FIG. 2: Theoretical enthalpy of formation vs. lattice constant 
for II-IV- V2 chalcopyrites. Negative enthalpy implies a stable 
host compound. 
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FIG. 3: Theoretical impurity formation energy of substitu- 
tional Mn vs. lattice constant of host chalcopyrite. Chemical 
potentials were defined by assuming the simultaneously Il-rich 
and IV-rich condition. 



that even substitutional Mn may be subject to a Jalin- 
Teller distortion, which can affect the magnetic interac- 
tions as discussed in Ref. pol |. 

It is important to realize that the impurity formation 
energy (and hence the solubility limit) in a multicompo- 
nent host is not uniquely defined, but depends on the 
growth environment — in particular, on the energy fi re- 
quired to transfer one atom from a reservoir into the 
host compoundSi. For the ternary chalcopyrites there 
are three such chemical potentials, {/in, /^iv, ^J'^/}^ which 
are related by + ^iv + 2/iv = Ai?/(II-IV-V2). The 
formation energies shown in Fig. 3 were evaluated for 
and /iiv at their highest allowed values, corresponding 
to a simultaneously Il-rich and IV-rich growth environ- 
ment. This condition can often be realized experimen- 
tally, unless a more stable compound intervenes. With 
this choice, half of the chalcopyrite hosts (including all 
of the carbides) prefer Mn on the group-II site, while 
half (including most of the silicides) prefer the group-IV 
site. For both the II and IV sites, there is the expected 
trend of larger lattice constants leading to smaller impu- 
rity formation energies — and hence to higher solubilities. 
This trend begins to disappear for lattice constants larger 
than ^^6 A; in this regime there is enough room available 
for the Mn impurity that the formation energies become 
essentially independent of both site and host. 

In the remainder of this work we assume that the 
prospects for ferromagnetism are best served when Mn 
substitutes on the group-IV site-'^^. The doping can be 
biased in this direction by changing the growth environ- 
ment to the simultaneously Il-rich and IV-poor condi- 
tion, which reduces the relative Mnjy formation energy, 
AE = Ai?/(Mniv) - AHf{Mnu), to its minimum phys- 
ically allowed value. The magnitude of the reduction de- 
pends on the nature of other possible intervening phases 
llxlVyVz, which we do not pursue here (see Ref. [231 for 



a related discussion). In the best-case scenario, the limit- 
ing IV-poor condition will lead to a change in AE by an 
amount AiJ/(II-IV-V2), relative to the values shown in 
Fig. 3. Using our calculated enthalpies, we find this lim- 
iting IV-poor condition favors Mn on the group-IV site 
in 50 of the 64 hosts (the carbides are again the domi- 
nant exception). Thus we conclude that under suitable 
growth conditions, substitutional Mn can be forced to 
preferentially occupy the group-IV site in the majority 
of chalcopyrite hosts. 

Even in such favorable cases, one expects some frac- 
tional occupation of Mn on group-II sites. Since this dop- 
ing is isoelectronic it is unlikely to produce either holes or 
electrons, and thus will probably require additional dop- 
ing to order ferromagnetically. A more pressing question 
is whether Mn doping on the group-IV sites necessarily 
leads to ferromagnetism in every II-IV-V2 chalcopyrite. 
And, for those hosts that can be made ferromagnetic, it 
is of great interest to know what Curie temperatures one 
can expect. 

These are very difficult questions to address compre- 
hensively, in a materials-specific framework, for any ran- 
dom magnetic alloy. Instead, we adopt again our earlier 
view and focus on a simple proxy for ferromagnetism, 
namely the nearest-neighbor Heisenberg coupling, J, for 
a specific distribution of Mn in each host. In the mean- 
field theory of Heisenberg ferromagnetism, the Curie 
temperature is simply proportional to |J|, with a pro- 
portionality constant depending on the spatial distribu- 
tion of Mn. Nevertheless, we caution that while qual- 
itative trends in J are likely to be meaningful within 
this materials class, the quantitative prediction of Curie 
temperatures in ferromagnetic semiconductors remains a 
difficult problem. Here we choose a fixed, experimentally 
reasonable Mn concentration of 12.5% (see, for exam- 
ple, Ref. [ni), and examine how J varies across different 
hosts. In a related study of the I-III-VI2 chalcopyrite 
CuGaSe2, Picozzi et al. showed that | J| increases with 
Mn concentration^^. Our method for computing J is the 
same as in Ref. "2^ , and gives nearly identical results for 
CuGaSe2 with 12.5% Mn. 

Surprisingly, we find that ferromagnetic alignment is 
favored (J is negative) in only half of the 64 chalcopy- 
rites. Three main factors determine the behavior of J 
across different hosts. The first is the group-II element: 
the majority of Mg- and Cd-containing hosts favor fer- 
romagnetic alignment, while the majority of Be- and 
Zn-containing hosts favor antiferromagnetic alignment. 
The second factor is the group-IV element: for many 
combinations of II and V elements, we find the order- 
ing J(C) < J(Si) < J(Ge) < J(Sn), except in the ni- 
trides (see Supplementary Information, Fig. SI). The 
group-V element plays the third — and lesser — role, with 
the ordering J(P) < J(As) < J(Sb) obtained for most 
cases except the carbides (see Supplementary Informa- 
tion, Fig. S2). None of these trends can be attributed 
to the approximate scaling, J ~ a~^, predicted by the 
mean-field solution of the Zener model24; indeed, across 
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this class of chalcopyrites we find no universal depen- 
dence of J on lattice constant (see Supplementary Infor- 
mation, Fig. S3). Nor do we find any simple correlation 
between J and the bandgap; such a correlation might 
have been expected if the small-gap hosts reflected Zener 
physics (with modest Curie temperatures) and the large- 
gap hosts double-exchange physics (with higher Curie 
temperatures)^^. 

Fig. 4 shows the computed values of J plotted against 
the relative Mnjv impurity formation energy, AE. We 
do not imply any correlation between these two quan- 
tities, but rather display them together because each 
is centrally important for designing new ferromagnetic 
chalcopyrites. In particular, the lower-left quadrant con- 
tains chalcopyrites for which occupation by, and ferro- 
magnetic alignment of, Mn on the group-IV site is fa- 
vorable under experimentally plausible II- and IV-rich 
growth conditions. This conservatively defined sweet 
spot contains just 19 chalcopyrites out of the full set of 
64. Of these, only eight are thermodynamically stable: 
four are lattice-matched to a mainstream semiconduc- 
tor (BeSnNa, MgGePa, MgSiAsa, MgCeAsa) and four 
are not (BeGeNa, MgSiSba, MgCeNa, MgCeSba). Seven 
of these eight have values of J comparable to, or larger 
than, that of CdGeP2 and thus can be expected to have 



similar Curie temperatures — of order room temperature 
or larger. This small set of materials should provide a 
propitious starting point for more detailed experimental 
and theoretical scrutiny. 

One can also consider the best-case scenario described 
earlier and consider AE in the limiting IV-poor growth 
condition. This leads to horizontal shifts, by ^Hf, of 
the points in Fig. 4 — leftward for negative enthalpies and 
rightward for positive — and thereby moves 13 more chal- 
copyrites into the lower left sweet spot. Ten of these are 
lattice- matched to mainstream semiconductors (MgSiN2, 
MgSiPa, MgSnPa, MgSnAsa, ZnSiPa, ZnGeNj, CdSiPa, 
CdSiAs2, CdGeAs2, CdSnP2), and three are not (BeSiP2, 
CdSiN2, CdGeP2); all are thermodynamically stable. 

Experimentally, CdGeP2 and ZnSnAs2 have ferromag- 
netic ground states with comparable Curie temperatures, 
while ZnGeP2 has an antiferromagnetic ground state (fer- 
romagnetism disappears below 47 K 11]). Our results, 
which presume Mn doping on the group-IV site, predict 
CdGeP2 to be ferromagnetic but ZnSnAs2 and ZnGcP2 
to be antiferromagnetic. Reconciling these findings is 
problematic, because the experimental group-II and -IV 
chemical potentials — and thereby the favored Mn site — 
are in general not known. The available evidence points 
to Mniv substitution in CdGeP2 [III, Mnn in ZnGeP2 
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[Tll |. and mixed occupation in ZnSnAs2 ^3- Future 
work addressing the magnetic interactions between Mn 
on group-II sites is in progress. 

Finally, we address the prospects for using ferromag- 
netic chalcopyrites as sources of spin-polarized current 
for injection into nonmagnetic semiconductors. Many is- 
sues play a role in this process, including the spin polar- 
ization of states at the Fermi level, the band velocities 
of those states, native defects and other structural inho- 
mogeneities, and spin scattering at the interface^. Here 
we focus on the first of these. For dilute Mn doping 
(12.5% and 6.25%, and presumably below this as well), 
all of the 21 stable ferromagnetic chalcopyrites identified 
above are "idealized half-metals," by which we mean at 
that at zero temperature and for an ordered arrangement 
of Mn, there are states at the Fermi level only in one 
spin channel. The effects of both finite temperature and 
site disorder will generally reduce the spin polarization 
from this idealized case. At higher Mn concentration 
(50%), half-metallicity persists for all but two of these 
21 hosts (BeSnNa, CdSiNa). Even when Mn completely 
occupies the group-IV sublattice, the resulting stoichio- 
metric compounds are half-metallic in all cases except 
for the nitrides. Thus the prospects for using ferromag- 
netic chalcopyrites as spin-polarized sources appear quite 
favorable. 

In summary, we have examined theoretically the 
prospects for ferromagnetism within the class of all pos- 
sible II-IV-V2 chalcopyrites having constituents chosen 
from the first four rows of the Periodic Table. As ex- 
pected for semiconductors, we find that the electronic 
properties of the host materials (band gap and enthalpy 
of formation) are closely related to their structural prop- 
erties (lattice constant, a). Contrary to conventional wis- 
dom, we do not find that Mn is generally more soluble on 
the group-II site than on the group-IV site; instead, the 
relative site solubility depends on both the nature of the 
host (e.g. carbide vs. silicide) and the choice of growth 
condition, while the absolute Mn solubility is largely de- 
termined by the host lattice constant. Our results for 
the Heisenberg coupling between Mn show that the Curie 
temperatures expected in Mn-doped chalcopyrites do not 
follow the scaling with lattice constant predicted by the 
Zener model; nor do they show any systematic varia- 
tion with band gap, as the double-exchange model would 
predict. By identifying those chalcopyrites that simul- 
taneously exhibit thermodynamic stability, favorable Mn 
doping site, and ferromagnetic Mn interactions, we iden- 



tify two small sets of chalcopyrites that show excellent 
prospects for stable ferromagnetism under realistic and 
attainable experimental conditions. 
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SpinS program. I.Z. acknowledges financial support 
from the National Research Council. Discussions with 
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Methods. All the numerical results reported here are 
based on density-functional theory in the generalized-gradient 
approximation^^, using ultrasoft pseudopotentials as imple- 
mented in VASP— ii^. For each host chalcopyrite, the plane- 
wave cutoff was separately determined by the three con- 
stituent elements (four if Mn was included) . For the undoped 
hosts, equilibrium lattice parameters and the internal coordi- 
nate were optimized using a 4x4x4 Monkhorst-Pack sampling 
of the Brillouin zone. Band gaps were evaluated by identify- 
ing band edges based on a 12x12x12 sampling that included 
the zone center and high-symmetry points. 

For enthalpies of formation, the total energies of the ele- 
mental phases were calculated using the following structures, 
which are either ground-state phases or energetically very 
close thereto: hep (Be, Mg, Cd, Zn); diamond (C, Si, Ge); 
white-tin (Sn); molecular dimer (N); black phosphorous (P); 
aAs (As, Sb). 

Mn impurity formation energies were computed using sim- 
ple cubic (or, for non- ideal c/a, tetragonal) supercells of 64 
host atoms. To properly represent the dilute impurity limit, 
all atomic positions were relaxed with fixed ideal supercell 
lattice parameters, within the constraint that Mn remain on- 
center. A single Monkhorst-Pack fc-point was used. Chemical 
potentials were defined with respect to their thermodynamic 
upper limit determined by the elemental phases listed above. 
For the Mn chemical potential, a non-magnetic fee structure 
was used to approximate the more complicated aMn ground- 
state phase. 

Heisenberg spin couplings for 12.5% Mn-doped chalcopy- 
rites were computed using the method described in Ref. |2^ 
The lattice parameters and internal coordinate of each su- 
percell were assumed to depend linearly on the Mn content 
(Vegard's law), and were interpolated using lattice parame- 
ters and internal coordinate calculated for the fully Muiy- 
substituted host. Two Mn per 64-atom supercell were then 
arranged on an uniform lattice and full atomic relaxation was 
performed. Total energies for two spin configurations were 
computed — ferromagnetic, and antiferromagnetic with [100] 
ordering wavevector — whereupon the energy difference gives 
J within the nearest-neighbor Heisenberg model. 
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Supplement 




C Si Ge Sn C Si Ge Sn C Si Ge Sn G Si Ge Sn 



FIG. 5: [FIG. SI] Theoretical Mniv spin coupling vs. group-IV element, for each of 16 sets of II-IV-V2 chalcopyrites. Results 
are plotted with arbitrary units to emphasize trends; the dotted lines denote J = 0. For many combinations of II and V 
elements, the ordering J(C) < J(Si) < J(Ge) < J(Sn) is found, except in the nitrides. 




FIG. 6: [FIG. S2] Theoretical Mniv spin coupling vs. group-V element, for each of 16 sets of II-IV-V2 chalcopyrites. Results 
are plotted with arbitrary units to emphasize trends; the dotted lines denote J = 0. The ordering J(P) < J(As) < J(Sb) is 
obtained for most cases, except in the carbides. 
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FIG. 7: [FIG. S3] Theoretical Mmv spin coupling vs. lattice constant of host chalcopyrite. No evidence is observed for the 
approximate scaling, J ~ a~^, predicted by the mean- field solution of the Zener model. 



